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We have studied the thermodynamics of isolated ho- 
mopolymer chains of varying stiffness using a lattice model. 
A complex phase behaviour is found; phases include chain- 
folded 'crystalline' structures, the disordered globule and the 
coil. It is found, in agreement with recent theoretical calcu- 
lations, that the temperature at which the solid-globule tran- 
sition occurs increases with chain stiffness, whilst the ^-point 
has only a weak dependence on stiffness. Therefore, for sufh- 
ciently stiff chains there is no globular phase and the polymer 
passes directly from the solid to the coil. This effect is anal- 
ogous to the disappearance of the liquid phase observed for 
simple atomic systems as the range of the potential is de- 
creased. 
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I. INTRODUCTION 

Most of the work on the thermodynamics of isolated 
homopolymers has concentrated on the collapse transi- 
tion between the high temperature coil and the lower 
temperature dense globule. This emphasis is a result of 
the extensive theoretical work on the transitionoH and 
the relative ease of simulations at these low densities — 
indeed simulations have been performed up to very large 
sizesj3 thus providing detailed tests of theoretical predic- 
tions. 

In contrast, there has been much less interest in the 
possibility of a low temperature order-disorder transition 
between two dense phases both because of the greater dif- 
ficulty of simulating dense polymers and because fewer 
theoretical expectations for such a transition are avail- 
able. This situation contrasts with that for heteropoly- 
mers where, in the context of protein folding, there has 
been intense interest in the transition between the molten 
globule and the native state of the protein. The structure 
of the native state reflects the amino acid sequence and 
the specific interactions between these units, and so it 
might be thought that an ordered structure is less likely 
when all polymer units are identical. However, this is 
not the lesson from other finite systems. For example, 
homogeneous atomic clusters show a rich low temper- 
ature phase behaviour; there is thejifijiite-size analogue 
of the first-order melting transition,aj3 surface melting,El 
and even low temperature transitions between different 
ordered forms .0 

Recently, the existence of an isolated homopolyjmp: 
order-disorder transition has begun to be confirmed.ETtHl 



In their study of a lattice homopolymer model which in- 
volved three-body forceS|-Kuznetsov et al. found phases 
with orientational order,Ll and in their simple oif-lattice 
model Zhou et al. observed a jOrder-disorder transition 
and also a solid-solid transitionJid Moreover, in some ear- 
lier studies of the collapsed polymer glimpses of these 
transitions were seen.ElTo Here, we add to this growing 
understanding of the low temperature phase behaviour 
of isolated homopolymers by studying a lattice model of 
a semi-flexible polymer, looking particularly at the effect 
of stiffness on the order-disorder transition. We com- 
pare our results jsdth the phase diagram calculated in 
recent theoretically and simulationE^I studies of the poly- 
mer model we use here. 



II. METHODS 

A. Polymer Model 

In our model the polymer is represented by an A^-unit 
self-avoiding walk on a simple cubic lattice. There is an 
attractive energy, e, between non-bonded polymer units 
on adjacent lattice sites and an energetic penalty, eg, for 



kinks in the chain. The total energy is given by 



E = -n 



(1) 



where Upp is the number of polymer-polymer contacts 
and Hg is the number of kinks or 'gauche bonds' in the 
chain, e can be considered to be an effective interaction 
representing the combined effects of polymer-polymer, 
polymer-solvent and solvent-solvent interactions, and so 
our model is a simplified representation of a semi-flexible 
polymer in solution. The behaviour of the polymer is con- 
trolled by the ratio kT/e; large values can be considered 
as either high temperature or good solvent conditions, 
and low values as low temperature or bad solvent condi- 
tions. The parameter eg defines the stiffness of the chain. 
The polymer chain is flexible at Cg^O and becomes stiffer 
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In this study we only consider eg > 0. 
eg = 0, this model corresponds to one originated 
, been much used to study homopolymer 
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The system with positive eg h. 



fias-been 

recently studied theoretically by Doniach et |flitJ and 
using simulation by BastoUa and Grassberger.E3 In our 
work we pay special attention to the structural changes 
of a polymer of a specific size. In this sense our work is 
complementary to that of ref. ^ which focussed on the 
accurate mapping of the phase diagram. 
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Our model was chosen because we wished to have the 
simplest model in which we could understand the effects 
of stiffness, and because it gives low energy states that 
have chain folds resptpbling those found in lamellar ho- 
mopolymer crystalsH Similar structures have previously 
been seerijn_a diamond-lattice model of semi-flexible 
polymers Jlirll3 and in moleculax jiynamics simulations 
of isolated polyethylene chains.t3~Cj Particular mention 
should be made of the studies by pK-olinski et al. on the 
effect of chain stiffness on coUapseEjiij since they found 
evidence for many of the phenomena that we explore sys- 
tematically here. 

The global minimum at a particular (positive) is 
determined by a balance between maximizing Upp and 
minimizing Ug. If the polymer is able to form a structure 
that is a cuboid with dimensions a x 6 x c (iV = a6c), 
where a < b < c, then 



(a) 



and 



= 2N - ab- ac- bc+1 



n"^"' = 2ab-2 



(2) 



(3) 



n™" have the 



The structures that correspond to Ug 
polymer chain folded back and forth along the longest 
dimension of the cuboid. By minimizing the resulting 
expression for the energy one finds that the lowest energy 
polymer configuration should have 



c 2e 

a = b and - = 1 H -. 

a e 



(4) 



Therefore, at = the ideal shape is a cube and for 
positive £g a cuboid extended in one direction, the aspect 
ratio of which increases as the chain becomes stiffer. As 
the ideal aspect ratio of the crystallite is independent of 
iV, its squared radius of gyration, Rg, will scale as N^^^; 
this scaling is the same as for the disordered collapsed 
globule. 

Substitution of the optimal dimensions of the cuboid 
[Equation (|)] into Equations gives a lower bound 

to the energy of the global minimum, 

Eopt/e = -2N + 37V2/3(1 + 2eg/e)i/3 _ i _ 2e<,/e. (5) 

However, at most sizes and values of eg it is not possible 
to form a cuboid with the optimal dimensions, and so the 
energy of the global minimum will be higher than given 
by the above expression. Nevertheless, it is easy to find 
the global minimum just by considering the structures 
which most closely approximate this ideal shape. 



B. Simulation Techniques 

Recent advances in simulation techniques have made 
it possible to begin to study dense polymer systenis. In 
particular, we use configurational-bias Monte Carlocj in- 
cluding moves in which a mid-section of the chain is 




FIG. 1. Bond-flipping moves, (a) Four-bond flip in which 
the mid-section of the chain is reordered. Moves of this type 
were only attempted if they preserved the integrity of the 
chain, (b) Two-bond flip which results in a new chain end. 
Only the bonds which change are explicitly depicted; other 
sections of the polymer are represented by wiggly lines. 



regrown. We also make occasional bond-flipping moves 
[Figure |l| which, although they do not change the shape 
of the volume occupied by the polymepnehange the path 
of the polymer through that volume.EIrEa These moves 
speed up equilibration in the dense phases. The simula- 
tion method was tested by comparme-tp results obtained 
for Cg — by exact enumeration. cSt3 Thermodynamic 
properties, such as the heat capacity, were calculated 
from the energy distxibiitions of each run using the multi- 
histogram method JEI'l3 

The number of Monte Carlo steps used in our simu- 
lations was typically 4-30x10^. The longer simulations 
were required for the larger polymers, especially at tem- 
peratures where two or more states coexisted. The sizes 
of the polymers we studied were N=27, 64, 100, 216, 343. 
These sizes are 'magic numbers' for eg = 0, because com- 
pact cuboids of low aspect ratios can be formed at these 
sizes. In the presentation of our results we concentrate on 
polymers with A^=100 and 343, the former because the 
smaller size allows for clear visualization of the structures 
of the different phases, and the latter because the effects 
of finite size will be smallest. 

In order to monitor the orientational order within the 
polymer we devised an order parameter, Q, which is given 

by 
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where is the number of bonds in the direction a and 
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TABLE I. Properties of the global optima for a 100-unit 
polymer, tg""' and e^""" give the range of €g for which a 
structure is the global minimum. c/\/ab is the aspect ratio 
and tg^^ is the value of the stiffness for which a crystallite 
with that aspect ratio is expected to be lowest in energy and 
is given by e^' = {c/Vab - l)e/2. 
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A 


136 


38 


0.000 


0.375 


1.118 


0.059 


B 


133 


30 


0.375 


0.500 


1.563 


0.281 


C 


129 


22 


0.500 


0.833 


2.406 


0.703 


D 


124 


16 


0.833 


2.167 


3.704 


1.352 


E 


111 


10 


2.167 


3.500 


6.804 


2.902 


F 


97 


6 


3.500 


12.000 


12.500 


5.750 





FIG. 3. Visualizations of various states of a 100-unit poly- 
mer, (a) A polymer based on structure D with some disorder 
in the stem lengths from a simulation at T = 0.9efc~^ and 
eg = 2e. (b) A folded structure with 8 stems that contributes 
to the middle peak of Figure ||b (T = 1.2efc~\ = 3e). (c) 
A typical configuration of the dense globule with i?g = 5.60 
(d) A typical configuration of the 
= 5.0efe-\ ta = e) 



(T = 0.75efc-\ eg = e). 
coil with Rl = 38.91 (T : 




FIG. 2. Global minima of the 100-unit polymer at different 
values of eg . The labels correspond to those in Table |. 



(TV — l)/3 is the expected value of Ua if the bonds are 
oriented isotropically. Q has a value of 1 if all the bonds 
are in the same direction, i.e. the polymer has a linear 
configuration, and a value of if the bonds are oriented 
isotropically. 



III. RESULTS 

A. Solid phase 

In Table | the properties of the global minima for 
A'^=100 are given for the range of eg we consider in this 
study, and examples of these minima are illustrated in 
Figure |[ The shapes of the global minima agree wel l 
with that expected based on the analysis of section [I A . 
In the final column of Table || we have given the value of 
the stiffness for which a crystallite with the aspect ratio 
of the global minimum is expected to be lowest in energy 
based on Equation (|^). This value generally lies in the 
middle of the range for which the structure is the global 
minimum. Most of the global minima have a degeneracy 
associated with the different possible ways of folding the 
chain back and forth. This degeneracy decreases as the 
chains become stiffer and the crystallites become more 
extended. For eg — there is no constraint on Ug, and so 



the global minimum has a much larger degeneracy, and 
the majority of the isomers of the global minimum have 
no orientational order. 

The global minimum is the free energy global minimum 
at zero temperature. However, as the temperature is 
increased it will become favourable to introduce defects 
into the structures. One of the mechanisms we commonly 
observed was through fluctuations in the lengths of the 
folds, moving in and out like the slide of a trombone. An 
example of such a structure is shown in Figure The 
generation of this type of defect is especially common at 
larger values of eg since it does not involve an increase of 

Ug. 

The dependency of the degeneracy, and therefore the 
entropy, on the aspect ratio of the crystallites raises the 
possibility of transitions to crystallites with a smaller as- 
pect ratio as the temperature is increased. Indeed, this 
is usually what is seen and leads to the decrease of Rg 
with temperature that is observed before melting [Figure 
^ . At the centre of the transition, where the free energy 
of each crystallite is equal, ideally, the polymer would 
be seen to oscillate between the two forms [Figure |a] 
during the simulation spending equal time in each. This 
'dynamic coexistence' of structures leads to a bimodal 
(or even multimodal if more than two forms are stable) 
distribution of Rg [Figure 0b] . Since the Landau free en- 
ergy is given by Ai(g) = A — kT\ogpq{q), where A is the 
Helmholtz free energy SLiiidpq{q) the canonical probability 
distribution for an order parameter q, the multimodality 
in Rg implies that there are free energy barriers between 
the different crystallites. In the example shown in Fig- 
ure ^, the peak with the highest value of R^ corresponds 
to structures similar to the global minimum D, and the 
peak with the lowest value of Rg have structures simi- 
lar to C. The other peak consists of structures that have 
eight aligned stems, either in a 4 x 2 array or a square 
array with one corner unoccupied [Figure ^] . 

However, isomerization between the different forms of- 
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FIG. 4. Behaviour of Rg as a function of temperature and 
eg for (a) iV=100 and (b) iV=343. Each hne is labelled by 
the value of eg/e. 
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FIG. 5. Solid-solid coexistence observed for a 100-unit 
polymer with eg = Se at T = 1.2efe~^. (a) Fluctuations in 
Rg during a 10 million step Monte Carlo run. (b) Probability 
distribution of i??. 



ten requires a large-scale change in structure, which is 
particularly difficult if these transitions occur at low tem- 
perature, where the free energy barriers to transitions are 
largest. With simulation techniques that only apply local 
moves, transitions between crystallites would be effec- 
tively impossible to observe and even with a technique 
such as configurational-bias Monte Carlo, where global 
moves are possible, transitions may be rare, particularly 
for the larger polymers in this study. This possible lack 
of ergodicity on the simulation time scales can lead occa- 
sionally to the abrupt jumps seen in (e.g. for A^=100 

and for A^=343 at eg=2e and 



at eg=4e and T=lek~ 

r=0.95efc~^) rather than the smoother transition that 
would be expected if equilibrium values of Rg were ob- 
tained. 

The differences in energy and entropy between crys- 
tallites of different aspect ratio are due to surface effects 
and so scale less than linearly with size. Therefore, these 
solid-solid transitions are not finite-size analogues of bulk 
first-order phase transitions. 

Interestingly, this coexistence of polymers with dif- 
ferent cuboidal shapes but the same basic structure 
bears some resemblance to the coexistence of cuboidal 
sodium chloride clusters that has recently been ob- 
served experimentallyBj All the clusters have the rock- 
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FIG. 6. Behaviour of (a) Q and (b) C„ as a function of 
temperature and eg for Af=343. Each line is labelled by the 
value of Eg/e. 



salt structure but the cuboids have different dimensions. 



B. Solid-globule/coil transition 

For Eg > the global minimum is oricntationally or- 
dered with a value of Q near to unity. Of course, this 
folded state corresponds to the free energy global mini- 
mum at zero temperature. Hovirever, as the temperature 
is increased there must come a point when a disordered 
higher entropy state, be it the globule or coil, becomes 
lowest in free energy, and the polymer loses its orienta- 
tional order — it 'melts'. This transition is signalled by 
a decrease in Q to a value close to zero, giving rise to a 
typical sigmoidal shape for the temperature dependence 
of Q [Figure This transition is accompanied by a 
peak in the heat capacity [Figure and we use the 
position of this maximum to define the melting temper- 
ature of the polymer, T™. (Alternative definitions, such 
as the temperature at which Q = 0.5, give practically 
identical results). The transition also often involves a 



change in the radius of gyration, the sign and magnitude 
of this change depending on the stiffness of the polymer. 
For lower values of Cg the transition to the dense globule, 
leads to a decrease in Rg, but at higher values of eg the 
transition to the coil leads to an increase in Rg {eg—Qt, 
8e in Figure |4^). 

At temperatures in the transition region, as for the 
solid-solid transitions, dynamic coexistence of the or- 
dered and disordered states is seen, leading to multi- 
modal probability distributions (and free energy barri- 
ers) [Figure |^ . In the example shown in Figure ^ three 
states are seen. Structures based on the global minima 
B and C give rise to the two ordered states. The low 
value of Q (~ 0.5) for the state associated with struc- 
ture C is a reflection of the considerable disorder that 
can be present in the solid state at Tm- The peak with 
the lowest Q value is due to the disordered globule and 
an example of a polymer configuration that contributes 
to this peak is given in Figure 

The effects of stiffness on this order-disorder transi- 
tion can be seen from Figures ^ and |^. The temperature 
of the transition increases with stiffness, and at larger 
values of eg this increase is a little slower than linear. 
Accompanying this change is an increase in the height of 
the heat capacity peak, and therefore the latent heat of 
the transition. This increasing energy gap between the 
ordered and disordered states is a result of the increasing 
energetic penalty for the larger number of gauche bonds 
associated with the disordered state. As Tm = Ai^/AS*, 
the dependence of the energy gap, A£^, on the stiffness is 
one of the main causes of the increase in Tm- This effect 
is reinforced by the changes in the entropy difference be- 
tween ordered and disordered states, AS": as Cg increases 
the number of gauche bonds decreases, thus lowering the 
number of configurations contributing to the disordered 
state. 

This behaviour is similar to simple liquids where de- 
creasing the range of the potential leads to an increas- 
ing energy gap between solid and liquid because of the 
increasing energetic penalty for the disorder associated 
with the dispersion of nearest-neighbour distances in the 
liquid, thus playing an ipjoartant role in the destabiliza- 
tion of the liquid phase.oEj 

The respective roles of energy and entropy in the in- 
crease of Tm with stiffness can be investigated more quan- 
titatively by approximating the free energy of the dis- 
ordered state by that for the ideal coil. This is a rea- 
sonable approximation, since the ideal expression for ng, 
^ideai _ 47V/(exp(/3eg) +4), fits the simulation values 
for the coil, and for the globule, fairly closely. The free 
energy of an ideal coil Aideai is 



A.. 



deal 



-NkT log (1 + 4 exp(-/3eg)) . 



(7) 



Decomposing this expression into its energetic and en- 
tropic components showed that the increase in the en- 
ergy with stiffness is the main contributor to the change 
in the free energy of the disordered state with stiffness. 
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FIG. 7. (a) Solid-globule coexistence observed for a 
100-unit polymer with = e at T = 0.675efe"^ 

Two-dimensional probability distribution in Q and R^. The 
labels B and C refer to the global minima of Table Q on which 
the structure of the polymers contributing to the maxima are 
based, (b) Crystal-coil coexistence for a 343-unit polymer 
with Eg = e at r = 3.45efc^^. Two-dimensional free energy 
profile in Q and E. The contours occur at intervals of O.bkT 
above the global free energy minimum. The contours over 
15kT above the global minimum are not plotted. 
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FIG. 8. Phase diagram of the 343-unit polymer. The phase 
diagram is divided into regions by the values of T,n and Tg. 
The position of the heat capacity peak associated with col- 
lapse, Tc, has also been included. The solid lines with data 
points are simulation results, and the dotted line is fro m the 
simple theoretical calculation of Tm outlined in Section III B . 



except at the higher temperatures relevant to the 'melt- 
ing' transition for larger values of eg. 

Use of Equation (Q) also aUows us to calculate a value 
of Tm if we approximate the free energy of the solid by 
only its energetic component [Equation (||)]. This sim- 
ple calculation gives surprisingly good agreement with 
the simulation data [Figure ||], especially at larger eg. 
It breaks down at low eg, e.g. the prediction of a non- 
zero value of Tm at eg=0, because Equation (Q) is the 
free energy of a coil in the absence of any interactions; 
for the dense polymers at low eg and near to Tm, the 
polymer-polymer contacts make a significant contribu- 
tion to the energy of the disordered state. The number 
of polymer-polymer contacts in the disordered polymer 
decreases significantly with increasing temperature, be- 
cause of the larger entropy of less dense configurations. 
Therefore, Equation (0) becomes a better approximation 
to the free energy at higher temperatures, and thus pro- 
vides a better description of melting at the higher tem- 
peratures relevant for larger eg. 

Our simulation results for Tm are also in qualitative 
agreement with theoretical results in—which a more so- 
phisticated treatment of the globuletj gives the correct 
behaviour for T,„ at low eg. However, a drawback of the 
description of Doniach et al. is that the free energy of 
the coil per polymer unit is nonzero in the limit of large 
eg/kT , causing Tm to reach an asymptotic value, rather 
than continuing to increase with eg. 

The increase in /S.E with stiffness also has an effect 
on the coexistence of ordered and disordered states. At 
large values of eg there is multimodality in the proba- 
bility distribution for the energy as well as for the order 
parameter, because of larger free energy barriers between 
the states. In the example shown in Figure ^ there is a 
free energy barrier of 3.24fcr for passing from the crystal 



6 



(a) 





\ 

\ 

\ 343 






\\ V .. 216 
















1 


; 27\ \ 




If 






Is 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
temperature / Ek ^ 



(b) 



343 




216 


ji 


100 n 


) : 1 


64 M 




TiriJij 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
temperature / eA' 

FIG. 9. Cv/N as a function of size for (a) = and (b) 



to the coil. 

A consideration of the effect of size on the transition 
shows that as the polymer becomes longer the melt- 
ing point becomes higher, and the transition becomes 
sharper with an increasing latent heat per monomer [Fig- 
ure |9|b] . This is consistent with the transition being the 
finite-size analogue of a first-order phase transition. T,„ 
increases with size because the effect of the surface term 
(the N^/"^ term) in the energy of a crystallite [Equation 
(^] diminishes with size. Moreover, as the coefficient of 
the surface term increases with eg (the higher aspect ra- 
tio crystallites have a larger surface area) the effects of 
size on are more pronounced at larger eg. 

At eg=0 the behaviour is qualitatively different be- 
cause there is no energy difference between the oricn- 
tationally ordered and disordered forms. As there are 
far fewer states that possess orientational order, they 
are never thermodynamically favoured and so there is 
no orientational order-disorder transition and Q always 
has a low value [Figure Despite this, the polymers we 
studied do have a low temperature heat capacity peak at 
T ~ OAek~^ [Figure Ba]. This feature stems from a tran- 



sition between the maximally compact cuboidal global 
minimum and a more spherical dense globule, and gives 
rise to bimodality in the canonical probability distribu- 
tion of the energy. However, the latent heat per atom for 
the transition decreases with increasing size, indicating 
that the transition is not tending to a first-order phase 
transition in the bulk limit. Furthermore, the form of the 
heat capacity can be very different for sizes at which it 
is not possible to form complete cuboids. 

The order-disorder transition obseipved by Zhou et al. 
was also for a fully flexible polymer.ll3 However, as the 
model used was off-lattice, unlike in our model, there is 
the possibility of a first-order transition associated with 
the condensation of the monomers onto a lattice. 



C. The globule to coil transition 

For polymers, there can be two disordered phases, the 
dense globule [Figure |^] and the coil [Figure ^jd] . They 
are differentiated by the scaling behaviour of the size of 
the polymer in each phase. For the globule oc N^^-^ 



N^/^. Between these two states 
ti 

is the 6'-point where Rg (x N and the polymer is said to 
behave ideally. On passing from the coil to the globule 
the polymer collapses, as seen in the sigmoidal shapes of 
the Rg curves [Figure [lO|a] . The collapse transition can 
also give rise to a high temperature peak in the heat ca- 
pacity, which is more rounded than that for melting and 
the latent heat of which is associated with the loss in 
polymer-polymer contacts on going to the lower density 



coil [Figure ||a]. These features are most clear at eg = 
where they are not obscured by the features due to the 
melting transition. At higher values of eg the sigmoidal 
shape of Rg is cut off at low temperatures by the rise 
associated with the transition to the crystalline states 
[Figure^, and the collapse transition just causes a high 
temperature shoulder in the heat capacity [Figure 
As the size of the polymer increases the transition be- 
comes sharper with a steeper rise in Rg [Figure [l0|a] and 
a narrower heat capacity peak at a higher temperature 
[Figure ^] . In the limit N ^ oo this heat capacity peak 
occurs at the 0-point. 

To complete the phase diagram for our model poly- 
mer we estimated Tg by making use of the scaling rela- 
tions: a series of plots of R'^/N should all cross at Tg. 
However, due to finite size effects the value obtained by 
this method only approaches the exact Tg (from below) 
as iV ^ cxD. For example, we estimate Tg{eg^0)~3.537 
from the crossing point of the lines for iV=216 andp343, 
whereas an accurate determination gives Te=3. 721.0 

Figure ^ shows that Tg has only a weak dependence on 
eg. Combined with the increase of with eg, this leads 
to the loss of the globular phase at eg ^ 12e. Above this 
value the polymer passes directly between the solid and 
the coil phases, i.e. there is only one disordered phase. 
An analogy can be made to the phase behaviour of sim- 
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FIG. 10. R^/N as a function of size for (a) €g — and (b) 



pie atomic systems, with a correspondence between the 
dense globule and the liquid and between the coil and 
the vapour. For both systems the denser phase disap- 
pears as an interaction parameter is varied, for polymers 
as the stiffness is increased and for atomics systems as 
the range of the potential is decreasedHj 

Our phase diagram is consistent with the work of Bas- 
toUa and Grassberger. From simulations of polymers 
longer than those we study here, theji, found that the 



globular phase disappears at eg ~ 13e.Ej The phase dia- 
grams is also very similar to the theoretical predictions 
of Doniach et al. who estimate the loss of the globular 
phase to occur at eg ~ 15eJl3 However, they assumed 
that Te is constant whereas both in ref. |l^ and here a 
small increase with eg is observed. 

Precursors of the loss of the globular phase can be seen 
at eg < 12. These effects are due to finite size, and for 
smaller polymers they occur at lower values of eg. Firstly, 
the features in the heat capacity due to collapse become 
engulfed by the peak due to melting. As the stiffness 
increases, the collapse peak becomes only a shoulder and 
then it disappears altogether [Figure Secondly, on 
melting the polymer can pass directly from the solid to a 
low density disordered state, albeit one where scales 
less than linearly with size. This effect can be seen for the 
27-unit polymer at eg = 3e [Figure |l^ , and is also seen 
for longer polymers at larger values of eg [Figure ^] . It 
results from the expansion of sufficiently stiff and small 
polymers with decreasing temperature, which is because 
the persistence length becomes a significant fraction of 
the total length. An estimate of the persistence length 
is given by the ratio N/ug, the average length between 
gauche bonds, N/nf""'- = 1 -I- exp(/3eg)/4. 

Considerable theoretical effort has gone into studying 
the evolution of the coil-globule transition as the stiffness 
of the polymer is increased, see Refs. ^ , p8| and references 
therein. The theories predict that the coil-globule tran- 
sition is second order foi^flexible polymers (this has been 
confirmed by simulationH) but becomes first order as the 
stiffness increases. Due to the small size of the polymers 
we have simulated we are unable to determine the order 
of the transition accurately and so cannot test this pre- 
diction. However, the theories all assume a liquid-like 
dense phase; they neglect the possibility of an ordered 
dense phase, despite it being well known that the 
state of DNA, a stiff polymer, has hexagonal orderEB 
It is clear from Figure || that for stiff polymers, the coil- 
globule transition is preempted by a coil-solid transition. 

The difference in phase behaviour between flexible and 
stiff polymers can be easily understood. As a flexible 
polymer is cooled, the Boltzmann weight for polymer- 
polymer contacts increases and so the number of con- 
tacts increases, at a cost in the entropy of the polymer. 
The entropy cost derives from the fact that the polymer 
must bend back on itself in order for the units to be in 
contact. The increase in the number of contacts is contin- 
uous and so the radius of gyration of the polymer varies 
continuously — the coil-globule transition is second or- 
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der. However, bending a stiff polymer back on itself is 
more difficult, it costs both entropy and energy. This 
larger cost can be repaid if more than one pair of units is 
in contact, which is true if the two parts of the polymer 
run parallel to each other for several units. Of course, 
as the polymer is stiff the entropy cost for two parts of 
the polymer to run parallel for a number of units is small. 
The low energy configurations of a stiff polymer are those 
with long parts of the polymer parallel, see Section HA. 
The energy gain is even l a rg er if a number of parts of the 
polymer form a bundle,cilc3 e.g. if four parts of the poly- 
mer which are running parallel form a square bundle the 
energy is not twice but four times that of two parts run- 
ning parallel. So, when the polymer is cooled below the 
point where the energy gain of bundles outweighs their 
entropy cost then these bundles proliferate and the radius 
of gyration drops suddenly — the coil-solid transition is 
first order. 

Our results can also be related to experiment. The 
coil-globule transition of polystyrene, a flexible polymer, 
is continuous and the globule is liquid-like.Ejo For DNA, 
an example of a stiff polymer, a first-order transition has 
been .observed between the coil and a compact dense 
state,E3a which has hexagonal order.cBc3 The lattice 
model used here also has a continuousEl coil-globule tran- 
sition when e^/e is small and a first-order transition from 
a coil to an ordered dense state when e^/e is large. It 
is thus able to reproduce the phenomenology of the ex- 
isting experimental data. However, we do not know of 
any examples where the three phases — coil, globule and 
crystalline — predicted by our model to occur for semi- 
flexible polymers have been experimentally observed. 



IV. CONCLUSION 

In this paper we have provided an example of an order- 
disorder transition that can be observed for an isolated 
homopolymer. Its basis in the interactions of the model 
is clear. The global minimum is orientationally ordered 
in order to minimize the energetic penalty for gauche 
bonds. At some finite temperature this order must be 
lost in a transition to the higher entropy globule or coil. 
As the stiffness of the polymer increases the energy gap 
between the ordered and disordered states increases con- 
tributing to an increase in the temperature at which this 
'melting' transition occurs. This effect, coupled with the 
weak dependence of the 0-point on the stiffness parame- 
ter, leads to the disappearance of the disordered globule 
for sufficiently stiff polymers; there is no longer a col- 
lapse transition from the coil to globule. This behaviour 
is analogous to the disappearance of the liquid phase of 
simple atomic systems as the range of the potential is de- 
creased. The phase diagram-ssje obtain is in g(>od agree- 
ment with recent theoreticalEj and simulationtS studies 
of this polymer model. 

One of the intriguing aspects of our model is the folded 



structures that form at low temperature. Firstly, these 
states are not artifacts of our lattice model. In simu- 
lations of isolated polyethylene chains, for the same en- 
ergetic reasons as in our model — the extra energy from 
polymer-polymer contacts outweighs the energetic cost 
of forming ..a . fold — relaxation to folded structures was 
observed. c3c3 Furthermore, as here, it was also found 
that the aspect ratio of the crystallites increased with in- 
creasing stiffness. However, we know of no experiments 
on a homopolymer system where evidence has been found 
for a transition from a disordered globule to a crystalline 
structure with chain folds. The most likely signature for 
such a transition would be an increase in the radius of 
gyration with decreasing temperature in very dilute poly- 
mer solutions. 

It is also natural to ask whether this study can pro- 
vide any insights into polymer crystallization, since poly- 
mer crystals have a lamellar morphology in which the 
polymer is folded back and forth in a similar manner to 
the folded structures we observe. However, the forma- 
tion of the folded structures for the isolated polymer is 
a purely thermodynamic effect — the structures are the 
globa.1 mipima — whereas in the bulk case it is a kinetic 
effectc3c3 — it is generally accepted that the global mini- 
mum is a crystal where all the chains are in an extended 
conformation. The study, though, may have a more in- 
direct relevance to the crystallization of polymers from 
solution. Although little consideration has been given to 
the structure of the polymer arriving at the surface in 
theories of polymer crystallization, it is not implausible 
that the existence of significant ordering in the adsorbing 
polymer would have a considerable effect on the crystal- 
lization process. 

Although the present study only examines the ho- 
mopolymer thermodynamics, some comments concern- 
ing the dynamics can be made. Klimov and Thiru- 
malai made the interesting observation that model pro- 
teins are most likely to be good folders when Tf/Te is 
large, where Tf is the temperature at which t^e transi- 
tion to the native state of the protein occursJlJ For our 
system T,„/Te increases with €g [Figure |^. Therefore, 
if Klimov and Thirumalai's relation also holds for our 
homopolymers, one would expect crystallization of the 
polymer to become more rapid as the stiffness increases — 
crystallization is easier direct from the coil than via the 
disordered dense globule. However, whereas at Tf there is 
a transition to a single state, the globule minimum, at our 
Tm a transition to an ensemble of crystalline structures 
occurs. Therefore, it does not necessarily follow that the 
homopolymers would be able to reach the global min- 
imum more rapidly as the stiffness increased — indeed it 
might be that trapping in low-lying crystalline structures 
is more pronounced for stiff polymers. 
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